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Phononless plasmon assisted transport through a long disordered array of finite length quantum 
wires is investigated analytically. Two temperature regimes, the low- and the high-temperature 
ones, with qualitatively different temperature dependencies of thermally activated resistance are 
identified. The characteristics of plasmon assisted and phonon assisted transport mechanisms are 
compared. Generically strong electron-electron interaction in quantum wires results in a qualitative 
change of the temperature dependence of thermally activated resistance in comparison to phonon 
assisted transport. At high temperatures, the thermally activated resistance is determined by the 
Luttinger liquid interaction parameter of the wires. 
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I. INTRODUCTION 

It is a well established fact that the single-electron transport through a one-dimensional system with random 
scattering is suppressed at zero temperature and infinitesimally small applied voltagei. The physical mechanism of 
the suppression of transport is known to be the Anderson localization phenomenon, which is explained as a constructive 
interference of initial and back-scattered waves that enhances the return probability after scattering by an impurity 
potential. Eventually the constructive interference leads to the localization of the single particle wave function in 
a one-dimensional disordered system. A natural conclusion follows, that suppression of the coherence of the single 
particle propagation leads to delocalization and favors transport through the system. Indeed, at finite temperature 
electrons couple to thermally activated bosonic excitations in the environment (usually phonons), which results in the 
decoherence of the electron motion and in the thermally activated electron transport!^. 

A basic theoretical question remains, whether the thermally activated transport can be induced by electron-electron 
interactions alone, that is without an external bath of bosonic excitations like phonons. The mechanism of a phononless 
thermally activated transport consists of dephasing of the electron wave function by electron-electron scattering that 
results from e-e interactions. In the case of decoherence by interactions, the role of bosonic bath is played by the 
electron-hole pairs, or charge fluctuations, that are excited thermally or in the result of electron-electron scattering. 
A microscopic mechanism for creating a bosonic bath out of localized electrons has been first proposed in Ref.^. 

Recently the problem of thermally activated transport in a disordered one-dimensional wire with interactions has 
been addressed in several paper8^i^*Si2iiS^. In the absence of disorder, a one dimensional wire with interactions is 
described by the Tomonaga-Luttinger liquid model. This model allows exact treatment of the low-energy physics in 
the interacting system using the bosonization technique^*^. At the same time it is known that the Luttinger liquid is 
destroyed virtually by any small amount of random scattering^. Therefore, in the localized regime that is dominated 
by disorder the Luttinger liquid description is inapplicable. 

The two abovementioned limits suggest two approaches to the problem of thermally activated transport. The first 
one is based on the description of one dimensional wire as a disordered Fermi liquid. In that approach, the interactions 
are treated perturbatively. In recent works a metal-insulator transition at finite temperature has been reported in a 
disordered weakly interacting systems. The mechanism of transport at finite temperature consists in formation of a 
many particle delocalized state in which correlations between distant localized single particle states are established 
through thermally excited multiple particle-hole pairs^iiSS. 

The other approach treats the localized system as a pinned charge density wave. The thermally activated transport 
is described as a propagation of instantons through the system^. The theoretical description has then much in 
common with bosonization. It is applicable for not too weak interactions. That description suggests a thermally 
activated behavior of conductance similar to the variable range hopping. Most importantly, an external bath of 
bosonic excitations with continuous spectrum is necessary to facilitate the transport^. 

The electron-electron scattering can be equivalently reformulated as a scattering of electrons by charge fluctuations 
(plasmons). If the effective interaction is short-range, the scattering rate by electron-plasmon collisions is given by the 
Fermi golden rule. The application of the Fermi golden rule is based on the assumptions of complete loss of coherence 
of plasmons between two subsequent collision events. This is a good approximation in a system with strongly screened 
interactions, which have a Fermi liquid ground state, such as a three dimensional metall. It becomes progressively 
worse with lowering the dimensionality of the system. An extreme example of a system with highly coherent plasmons 
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is given by a pure one-dimensional system, or quantum wire. In that case, the charge density fluctuations are well 
defined quasi-particles with an infinite coherence time, and theoretically such system is described by a Tomonaga- 
Luttinger (TL) liquid modetfi^. At the same time, the electron itself ceases to be a quasiparticle that propagates 
coherently through a Tomonaga-Luttinger liquid. By entering the TL liquid form outside, an electron decays into 
charge and spin excitations. This decay process has a very slow power law character in time, which is also known as 
a demonstration of the Anderson orthogonality catastrophei^i. For such systems a straightforward application of the 
Fermi golden rule leads to incorrect results as it will be shown on an example considered in the present paper. 

In this paper we investigate phononless thermally activated transport through a quasi-one-dimensional system 
formed by a parallel arrangement of conducting wires (see Fig. ^ . Each wire has a finite length L and the transport 
direction is perpendicular to the wires. We identify the regime, where the localization length of plasmons in the array 
exceeds very much the single particle localization length. In that regime we show that charge-density fiuctuations 
(plasmons) in the array can act as the agent promoting thermally activated transport, thus providing the possibility 
for phononless inelastic transport. The plasmon assisted transport mechanism considered in the present paper is 
essentially the same as in RefsiSliSS. However, the adopted model exhibits a crossover from a Fermi liquid metal (or 
Anderson insulator in the presence of disorder) at zero temperature to a strongly interacting sliding Luttinger liquid 
system at higher temperatures. Therefore, perturbative treatment of interactions at finite temperature is inapplicable 
to the system under consideration. As the result of generically strong plasmon-electron coupling in a quantum wire, the 
features of plasmon and phonon assisted transport are qualitatively different. We provide a qualitative explanation 
of plasmon assisted transport, identify the transport regimes, where the features of plasmon and phonon assisted 
transport are either similar or substantially different, and derive analytic expressions for the temperature dependence 
of the thermally activated resistance. 

Hopping transport in the considered model is of much relevance to a number of experimental setups, including 
quantum wire arrays in heterojunctions"'^^, carbon nanotube films-^^, atomic wires on silicon surfaca^^, and stripe 
phaseai^. At finite length of constituent wires, such systems represent particular examples of granular arrays, where 
a one-dimensional wire plays the role of a grain. Considered as a granular array, the array of parallel quantum wires 
is rather peculiar because of the very long charge relaxation time in a one-dimensional wire. Due to this peculiarity, 
the theoretical description of thermally activated transport in arrays of long quantum wires requires taking into 
account the charge dynamics and treatment of interactions beyond the capacitive model adopted in recent theoretical 
investigations of transport through disordered granular arrays ^i^^i^^ . 

The rest of the paper is organized as follows. In Sec. the theoretical model is introduced in detail and the 
relevant physical regime is identified. In Sec. Illll we derive general equations for the thermally activated resistance of 
the array. Results of detailed calculations in the low- and high-temperature regimes are presented in Sees. II VI and fVl 
respectively. In Sec. I VII the calculation of the thermally activated phonon assisted transport for the considered model 
is presented in order to make a comparison with the plasmon activated transport. Our findings are summarized in 



The model we formulate below is special, because it combines two seemingly incompatible features: i) it is strongly 
disordered for single electron transport; ii) it is much weaker disordered for propagation of plasmons. 

Consider a one dimensional array of parallel identical quantum wires of length L and diameter a placed regularly 
with the interwire distance d, L ^ d ^ a. We investigate transport in the direction perpendicular to the wires (see 
Fig. 1). The spectrum of low-energy plasmons in a single isolated wire is equidistant with energies 



where L is the length of the wire, and Ui is the plasmon velocity along the wire i. For identical wires at regular 
positions, the intra- and interwire interactions between the charge density fluctuations do not change along the array. 
Then each localized plasmon level broadens into a plasmon band with truly continuous spectrum, quite analogously 
to the formation of electronic bands in the tight binding model. The role of the hopping in the tight-binding model 
for plasmons is played by the matrix element of the charge-density interactions between the neighbor wires. The 
formation of plasmon bands is reflected by the dependence of the plasmon velocity along each wire on the wave vector 
p along the array (that is, perpendicular to the wires) Ui u{p). The particular form of plasmon dispersion depends 
on details of the interwire interactions, yet the function u(p) should be periodic with a period of one Brillouin zone. 
That is why we choose the specific form 
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II. THEORETICAL MODEL 
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Up = vq + vi cos(7rp), —1 < p < 1. 
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FIG. 1: Geometry of the model. The arrow shows the direction of the current. 



The chosen form of dispersion can be considered as the first two terms of the Fourier expansion of some general 
dispersion law. The plasmon energy within a band centered around the level n is given by en(p) = "t^'^- 



A. Single particle localization length 



Since we are interested in the single particle transport in the direction perpendicular to the wires, we consider in 
this section only the transverse components of single particle wave functions, assuming the factorization of the wave 
function into the transverse and longitudinal (with respect to the direction of the wires) parts. 

Suppose the typical height of potential barriers between the neighbor wires is very large. Let us describe the array 
without interactions by a tight-binding model. Then a single particle wave function has most of its weight inside the 
wires. We approximate a Vannier wave function localized in a single wire with number n by the form 



My)^^o{o{y+^- nd) + 9 (rid + ^ - v) + e-^l"'^-^l [e (nd - ^ ^ y) + 6 (y ^ ^ ^ nd)]] , 



(3) 



where 9{y) is a step function. Here and everywhere further y denotes the coordinate in the direction perpendicular to 
the wires, that is in the transport direction. Using the normalization condition for ipniy), the value of |V'oP is found 
to be 
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In the absence of disorder and interactions the single particle motion is described by very narrow energy bands with 
dispersion efe = tcos{kd), (— ^ < k < ^). With the help of ||2J) the bandwidth t that equals the hopping energy in the 
equivalent tight-binding model can be estimated as 
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where A describes the decay of the single particle wave function inside the barrier of the height H, 



A : 
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■H. 



(5) 



(6) 



Now let us introduce the disorder as a random height of the energy barrier between the neighbor wires. Furthermore, 
the charge transfer between the wires is accompanied by Coulomb blockade effects. Due to the random environment 
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(local concentration of charged impurities) around each wire, the charging energy Ec is random, which represents 
another source of randomness in the model. Such disorder induces fluctuations of the decay parameter A thus 
rendering t random. As the result, the single particle wave functions become localized. In ID the single particle 
localization length is of the order of the mean free path If ~ vpTf. The mean free time Tf is in turn related to the 
fluctuations of the heights of the tunneling barriers 

{SH{n)SH{n')) ^ ^ y 5„,„>. (7) 

ZTTiyiTfCL 

Here n and n' denote the numbers of the wires adjacent from the left to the barrier, and we assigned a length scale 
1/d to the (5-function. Furthermore, ly'i is a single particle density of states, which is given by j^i = in the center 
of the band. Now let us establish a connection between the mean free path / / and the fluctuations SX of the decay 
parameter of the wave function. Using the relation © we obtain 



6X = J^^SH. (8) 



Further we solve ITJ for Tf and express 5H through with the help of (jSj). In the result we obtain 

lf = 5 — = . (9) 

^ 8^((5A')(aA + l) ^ ' 

Imposing the condition for the single particle to be localized within two nearest neighbor wires lf = d,we obtain for 

B. Plasmon localization length 

The random height of interwire tunneling barriers affects the interwire charge density interactions through the 
randomness of the transverse part of the wave functions ipi. The matrix element of the charge density interaction 
between the neighbor wires with numbers 1 and 2 (direct part of Coulomb interaction) is calculated as 

, , \Myi)\'\My2)\' 

dyidy2- 

l"^- y2 -~ yi\ 

-a/2 /.a/2 e'^^l^^ I | ^2 (2/2 ) | ' 




a/2 J J -a/2 \d + y2 ~ Vl] 

yi\ \J-oo Ja/2 \d + y2-yi\ 



(11) 



We took into account that the wave function decays exponentially inside the barrier (for \yi\ > a/2). For Ad ^ 1 the 
main contribution to the matrix element is given by the first term in Hll|l that relates to the interaction of charge 
densities inside the wires. That term results in 

~a2|^,nv>,+i|Vrf, (12) 

where 
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IV'.r - - / dy\Uy)\- (13) 
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Note that being integrated over the wire crosssection, the value in insensitive to random variations of the wave 
function across the wire. The contributions from other terms in are suppressed as 1/A. The effect of the random 
potential barrier on the quantity is encoded in the decay factor A^. Setting in Q Ai = A + SXi we obtain the 
fluctuation of in the form 



(14) 
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Furthermore, using (|12|l we estimate the fluctuation of the matrix element of the interaction as 

SV..,,^^^. (15) 
In the bosonized form, the density-density interaction between the neighbor wires reads 

His+i ^ dxVi^i+i {dx(t)i{x)) {dx(l)i+i{x)) . (16) 



Using the mode expansion of the bosonic fields 4>i{x), 0^+1 (a; )i 

^,ix)^^Y.- (e'^'^k, + e'^^,) (17) 



together with the relations q = j^fig, cuq = uq, where u is the phase velocity of a plasmon along the wire, we rewrite 
(|16|l in the form of a hopping Hamiltonian for the bosons represented by operators bi^q, b\q 

q 

In the result of random interwire interactions, the plasmons become localized in the direction along the array. The 
plasmon localization length can be evaluated as = UgrTp , where denotes the mean free time of a plasmon mode 
with frequency uj, and Ugr is the group velocity of plasmons. The plasmon mean free time is related to the fluctuations 
of the matrix element of the interwire charge density interactions 

where z/p — ^a^^„ is the plasmon density of states in the middle of the band n. In turn, the velocity vi is related 
to the average value of the charge density interactions V — (Vi.i+i) that in turn characterizes the width of plasmon 
bands in the pure system. From the dispersion law of the lowest plasmon band = ^(wq + ^^i cos(7rp)) we infer 

V=^v,. (20) 

Furthermore, the plasmon group velocity is given by 

1 dtp TT^d 

Using H12|l and Q the average interaction strength V can be expressed as 

where ((5A^) ^ A'^ is assumed. Solving p9|l for and expressing the plasmon density of states for the first plasmon 
band n = 1 through the average interaction strength V with the help of (|20|l and (|22|l . we obtain the expression for 
the plasmon mean free time in the form 

u x2 rf(aA + l)4 
2nLjJ 8((5A2)a2 " ^ ' 

Now let us evaluate the ratio between the plasmon localization length Ip and the electron localization length d. 
The condition of strong electron localization on the length d determines the value of (SX^) as given by Eq. 1)1(1 |l . 
Substituting that value of (SX^) into the expression H23I) and multiplying by Ugr/d with Ugr given by (|21() we finally 
obtain 

d Alv^X d^ 
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This result agrees qualitatively with the evaluation of plasmon correlation length in a randomly inhomogeneous 
Luttinger liquid by Gramada and Raikh^^. Eq. H24|l shows that the condition Ip/d ^ 1 can be formulated as 



\d > In 



2uj\d 



s(aA + l)2 



(25) 



Eq. H25() defines the regime, where the plasmon localization length exceeds very much the electron localization length. 
If, additionally, Ip exceeds the length of the array in the transport direction (orthogonal to the wires), the plasmon 
bands can be applied to describe the plasmon spectrum in the array. 



III. RESISTANCE OF A DISORDERED ARRAY 



The resistance of a long disordered one-dimensional array is determined by so-called breaks, the junctions between 
two neighboring wires with exponentially high resistancei^. Let us denote the energy cost to transfer an electron 
over the break as Ea- To facilitate the transport over the break, the energy Ea should be borrowed by absorption 
of a bosonic excitation. For two isolated wires forming the break, the matching condition between the energy of a 
plasmon mode n in the wire i given by 1^ and the energy Ea, i?i,n = Ea, cannot be satisfied for arbitrary Ea because 
of discreteness of Ei^n- However, due to the charge-density interwire interaction (|18|l the energy can be transferred 
between excitations localized in different wires. In the regime discussed in previous section the plasmon localization 
length exceeds the length of the array. Then the plasmon band description of the plasmon spectrum is applicable. 
Treating the interwire charge density interaction perturbatively, we can write the transition rate with the absorption 
of a plasmon similarly to a transition with the absorption of a phonon using the Fermi golden rule^, 

7 (X j dpY, \Vn{p)\NB{en{p))f i-E^) [1-f {Ea + Ek)] S (e„(p) - Ea - E„, ~ E^) . (26) 

n m,fc— 

Here Vn{p) is the strength of interwire charge density interaction for the plasmon mode n, Nsiinip)) is the occupa- 
tion number of the plasmon mode, f{Em) denotes the Fermi distribution and describes the occupation of the m-th 
single-particle energy level in the wire, Em — ^^m. For narrow plasmon energy bands, the perturbative approach 
suggests that when the energy Ea lies in the gap between the plasmon bands the hopping over the break is blocked. 
This suggestion turns out to be wrong because of a conceptual difference between the plasmon and phonon trans- 
port mechanisms. Whereas the phonons represent a bath of bosonic excitations that is independent of electrons, the 
plasmons are "made" of electrons themselves. Consequently, while the electron-phonon interaction can generally be 
treated perturbatively, the perturbative treatment of plasmons is possible only under special conditions. The appli- 
cability of the perturbative treatment of plasmons is determined by the relation of two time scales: the characteristic 
time of plasmon dynamics tp and the characteristic time of a single electronic hop t^- If the Coulomb interaction 
in a grain is well-screened or the plasmons are strongly localized, then tp is the characteristic relaxation time of a 
plasmon excitation within a single grain. For tp <^ th the plasmons can be neglected in transport. The description 
of interactions thus reduces to the capacitive modelSiiSii^. For the delocalized undamped plasmons, the time tp is 
associated with the formation of an extended in space plasmonic excitation. In that case, tp <^ th correspond to the 
regime of a strongly nonlinear coupling between plasmons and electrons, and the perturbative treatment of plasmons 
is incorrect. Plasmons in one-dimensional wires represent a profound example for that regime. In particular, the 
relation tp <^ th is always fulfilled at the break. For the model considered in this paper tp ~ L/t;i. As we show below, 
due to the strong electron-plasmon coupling, the nonlinear effects lead to the creation of plasmon complexes with 
energies covering the whole spectrum continuously, even though the plasmon bands initially are very narrow. This in 
turn leads to plasmon assisted transport with a temperature dependence that is qualitatively different from the case 
of phonon assisted transport. In the regime tp ^ th, the effective interaction time is limited by th- Then the plasmon 
dynamics is essentially independent of the electron dynamics and does not affect the electronic transport. 

The resistance of the array is calculated along the lines of Ref>i&. Let us parameterize the tunneling matrix element 
between two wires in the form i^^i+i = exp(— jj/i^i+il/d). The parameter y can be associated with an effective distance 
between the two wires. This effective distance is random, its distribution follows from the distribution of the heights 
of potential barriers. Since a break, being a junction with exponentially large resistance, is not shorted by other 
resistances connected in parallel, we can write the resistance of a break in the form 

i?i =i?oexp[2|y,,,+i|/d + /(Sa,T)]. (27) 

Here Ea denotes an addition energy to transfer an electron over the break. We remind that the disorder enters the 
model only as a random distribution of addition energies Ea and effective distances yi^i+i- The function f{Ea,T) 



7 



accounts for the effect of thermally activated plasmons in the resistance of the break. According to Ref.J^, the 
probability density p{u) for the resistance R/Rq — e" is proportional to e~^^, where A is the area in the (y, Ea) phase 
space that results in the resistance e", and g is the linear density of localized one particle states. The area A can be 
calculated as 



Aiu) 



u-f{E,T)=0 



dE[u-fiE,T)]. 



(28) 



Then the resistance is given by 



R = Roly 



du 



(29) 



where ly is the length of the array. The function e"~^"^(") is usually strongly peaked, which justifies the evaluation 
of H29|l by the saddle point method. The saddle point value of u corresponds to the resistance of an optimal break 
in the chaiiii^. Therefore, in order to calculate the resistance of the array in the localized regime, we have to obtain 
an expression for the resistance of the break Ri. Since the break is not shorted by other resistances, we conclude 
R\ = l/tTi, where cti is the conductance of a break. Assume that the break is formed by a junction between the 
wires with numbers and 1. We take the position of a pinhole connecting the two wires as x. In the linear response 
approximation the current through the break Ji is determined by the correlation functionSS 



X{t) = (Tr (*o(a:,r)^'l(a;,r)^'i(x,0)*t(a;,0) 



(30) 



that characterizes the probability of a single hop over the break. Here ipi is the tunneling matrix element, r is an 
imaginary time. In what follows we denote X+{t) = X{t > 0) and X-{t) = X{t < 0). 

A unique feature of the chosen model is the applicability of the bosonized description that allows exact treatment 
of interactions and hence nonperturbative treatment of plasmons. Precisely, the plasmon dynamics is described by 
the action 



S 



dp 



dr 



L/2 



dx 



L/2 '^Kp 



1 



(31) 



representing a finite size generalization of the sliding Luttinger liquid modeUS. The relation of the plasmon velocity Up 
and the Luttinger liquid constant Kp with inter- and intrawire interactions has been calculated in Ref.^^. A fermion 
annihilation operator in the wire n, \['„(a;), is represented as 



i / dp(j)p{x)e 



(32) 



where x — R,L denotes the chirality, (f>p{x) is a chiral bosonic field, and is a Klein factor. The chiral field cj)^ 
in turn expressed through the field Qp{x) and its dual $p(a;). 



IS 



<^^(x) = (ep(x)±<&p(x))V^. 



(33) 



In the bosonized representation (|32|l . the correlation function X{t) factorizes in the correlator of Klein factors and 
the correlator of bosonic exponents that we denote as Xii{t). The time dependence of the Klein factors Fn.x{T) is 
given by the ground state energy of the wire n that includes the capacitive interaction between the wires. Thus the 
correlator of the Klein factors is proportional to e~^'='^, where Ec denotes the charging energy. Denoting the correlator 
of bosonic exponents as 



X^^ 



(r) = (t. (e 



-40Ox(O,O) 



)) 



we can cast the expression for the current into the form 



dt 



E 



,(^+B.)t^xx'(^ = i< + 0) - e*("-^=)*X»^'(T = it + 0) 



(34) 



(35) 



Here Xi x' = denote the chirality of corresponding boson field. We assume the size of the tunneling region along 
the wire much larger than the Fermi wave length. Then the terms with equal chiralities, x = x'j the major 
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contribution to the current 1)35(1 . Leaving only those terms and noticing that the contributions form the modes with 
left and right chiralities to the tunneling current are equal, we can write the conductance of a single junction as 



- ^^^^^^ 



(36) 



Substituting the explicit form of the correlator of free bosonic fields, assuming a symmetric form of the dispersion 
u{p) = u{—p), we finally cast the correlator (|34() to the form (0 < r < /?) 



Xb+{T) ^ exp [- {KpSp{T, a))] 



(37) 



where Kp = Kp + 1/ Kp, 



oo 

Sp (r, a) — I In sinh 



m=0 

lnsinh|^((m+l)/3-a) 



^((m+l)/3-r) 



In sinh 

and the average over the plasmon wave vector p is defined by 



TTM- 

"2L 



In sinh 
^(TO/? + a)l| 

J J J 



^(m/3 + r) 



)p= ■ (1 ^ 2cos(7rp) + cos(27rp)) dp. 
Jo 



(38) 



(39) 



For — /3 < T < the correlator assumes the form 

Xb- (r) = exp [- {KpSpi-T, -a))] . (40) 

In Eqs. 1)37(1 . 1(38(1 . 1(40(1 . the zero temperature result is given by the terms in 1(38(1 with m = 0, whereas the terms with 
771 > give the finite temperature correction. 

Further we assume the coupling constant to be p-independent, Kp = k and use the simplified dispersion law |(2Jl. 
Note, that in approximate evaluations it is much more important to keep the p-dependence of the velocity Up that 
reflects the formation of plasmon bands than the p-dependence of the coupling constant Kp. To the lowest order 
in p-dependent terms, the latter just leads to the averaging of the single Luttinger liquid result over the coupling 
constant. 

Expanding the logarithmic functions in ((38(1 we can rewrite 13711 . ((4011 in the form 



m— 1 



where 



(41) 



(42) 



Here and everywhere further "reg" denotes the short time regularization, which for a function /(r) is defined as 
follows 



Z„,(r) = exp i K ^ /e-^("'^+^) + e-^'^"'^-^ 



f{ry-^ = /(r) - /(a) 
a being a short-time cutoff. The function I± (t) is defined as 



tO / \ J T^{Up)p 



2L 



rog 
P 



(43) 



(44) 



IV. RESISTANCE AT LOW TEMPERATURES 



The low temperature transport regime is determined by the condition that the temperature is less than the distance 
between the neighbor plasmon bands, T < ttvo/L. In terms of 1(41(1 . ((42() and 1(44(1 the current through a break 1(35(1 
can be represented as a convolution in frequency space 



27re f°° 

- — |toi|' / duj [l"_{uj + Ec- eV) -ll{uj-Ec- eV)] Z{-uj), 



(45) 
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where denote the real time Fourier transform of the function and Z{u!) is the Fourier transform of 

■^(^) ~ 11™=! ^rn{t)- One can consider 1451) as a kind of Fermi golden rule formula by rewriting it as 

h{V) = e{T+{V)-r_{V)) (46) 

with obvious identifications for the transition rates 



r± = 



— ItoiP / dLulliLOTE,~eV)Z{-Lu). (47) 



Then the quantity Ij_ describes the tunneling density of states at zero temperature, while the factor Z{uj) describes 
the density of thermally activated plasmon complexes that play the role of bosonic agents assisting the hop. 

Expressions l|42|) . (|44|l contain the averaging over the plasmon wave vector p that mathematically describes the 
influence of plasmon bands on transport. After the analytical continuation to real times t = ir, the basic average to 



be used in subsequent calculations reads 



loibvi) + {1 - -L)i^{bvi) 
bvi 



(48) 



where b = (m/3 ± it)Tr/L, and Ii,{z) denotes the Bessel function of complex argument. The time t in H48|) is limited 
by the hopping time over the break t^. In what follows we concentrate on the regime of long hopping times th that is 
defined by the condition 1/th <C In that regime the charge density fluctuations from spatially separated regions 
can reach the break during the hopping time. Thus, the full plasmon spectrum can participate in the creation of a 
bosonic excitation that assists the hopping process. Mathematically, this fact is reflected in using the approximation 
of Bessel functions at large value of the argument in the form 



e 



y 7r^|wi|(TO/5 ± it) 



where w = vq — vi. Despite being obtained for vi < vq, 1)49(1 is essentially nonperturbative in vi , which reflects a highly 
nonlinear influence of plasmons on the thermally activated transport. We note in passing that the consideration of 
the long hopping time regime is in accord with the identification of the break as a junction with an exponentially high 
resistance. 



A. Zero temperature tunneling density of states 



Performing the analytical continuation to real times in H44II and using the large time approximation (|49|) for m = 0, 
we obtain ll{t) as an expansion 



2L 



ll \±nT^\vi\t 



1/2 



/ oo ^ 
\n—l 



(50) 



In performing the Fourier transformation in (|50|1 term by term, one meets the fohowing typical Fourier transform 

(51) 



dt- 



-i(i±io)]"/2 



2 sin 



2 yr^l--)0(Tc^)|u;|"/2-i 



It follows from H51|) that each term in the Fourier transform of H50|) is proportional to a step function of the type 
9 {(jj — Ea — 2i|h)^ the integer k being given a sum of I numbers of ni, ...,71; from the sum over n in jSUJ. Since 
Lo = eV, only the terms with lowest / and n are important in the linear transport regime at small uj. Leaving only the 
terms with ^ = 0, 1 and n = 1, we obtain 



I 2L 
7r2|vi| 



9 ±uj 



TT 



1j 



-1/2 



(52) 



where we introduced the addition energy Ea = Ec + ^{up)p- Appearance of a continuous tunneling density of states 
at finite bias uj — eV > Ea + ^''^ indicates a strongly correlated nature of the system under consideration. It is due 
to the highly nonlinear interaction of electrons and charge density waves. This feature, which is not relevant for zero 
temperature linear response because of the finite energy gap to the continuous part of the spectrum, will determine 
the thermally activated transport that is described below. For the zero temperature tunneling density of states we 
restrict ourselves with the zero approximation leaving only the delta- function term in H52|) . 
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B. Influence of thermally activated plasmons. 



As it has been pointed out above, the influence of thermally activated plasmons is contained in the factor Z{uj) 
31 that is given by the Fourier transform of 



Z{t) = exp 



fj— ±1 m— 1 



(53) 



At low temperatures, T ^ the major contribution to Z{t) is given by the term with m = 1. Leaving only that 
term, expanding the logarithms in (|53|l and using ()49|l we obtain 



Z{co) « £ 



n,i=0 



nW. 



2L 



-;3(^(n+0 + |--.,„t)|^_g^^J 



._isin(^)r(l-f) 



(2/3)- 



(54) 



where ein(w) = ^{l — n), and ly = l,n for Lo — ein ^ respectively. Each term in (|54|l describes a thermal excitation of a 
multiparticle plasmon complex with a continuous density of states. The plasmon complexes described mathematically 
by (|54|l form the bath of bosonic excitations that facilitate transport over the break. Higher values of m in (|53|l would 
correspond to excitations involving progressively more plasmon modes. The leading terms in (|54|l at low temperatures 
are given by ^, n = 0, 1. Those terms result in the leading low-tempcraturc contribution to cri in the form 



0-1 



-'LT (1) 
h \/27r2|ui|T 



'2\l3Ea\ 



(55) 



from which wc identify 



2n 



f{Ea,T)^P[\Ea\ + —w] -\n[-\l3Ea 



1, 



-3/2 



\f3Ea\ 



-1/2 



(56) 



C. Resistance of a long array 



Further evaluation of the resistance proceeds according to Eqs. (|27|) . (|28|) . (|29|) of Sec. Illll The upper limit of 
integration in l|28|l 



u-fiEo,T)^0 



(57) 



relates the value of u to the corresponding value of the addition energy Eq. Under the assumption (3Eq ^ \n{(3Eo) 
that will be selfconsistently confirmed by the solution below, the integration in H28|l results in 



(58) 











\2T 









The saddle point value of u in l|29() relates to the saddle point value of the energy Eq which is obtained as 



(59) 



Therefore, the assumption above is justified for low temperatures, when l/{Tgd) 3> 1. Further evaluation of 129() 
leads to the following result for the resistance of the array 



R 



h n^\vi\T^ 



■ exp 



1 



2gdT LTJ 



(60) 



For comparison, as it will be shown below in Sec. IVII in the case of phonon-assisted transport in the same model the 
preexponential factor in H60|l goes like T^/^. 
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D. Low temperatures, short array. 

In the case of a short array, the optimal break cannot be found in the array. In that case the highest resistance 
amongst all those present in the array determines the total resistance. That highest resistance is defined by the 
condition lyp{uf) ~ 1, where ly denotes the length of the array, and p{u) = e~^^^^^ denotes the probability density 
for the logarithmic resistance u = log(i?/i?o)^- The value u/ is related to the length of the array by the condition 



1 



■ hiL, 



(61) 



Furthermore, the condition for the upper limit of integration Uf ~ f{Ef,T) in H28|l determines the characteristic 
value of the addition energy Ef. The area A{uf) is then given by H58I) with Eq replaced by Ef, that together with 
(|61|l provides the relation between the length of the array ly and the characteristic energy Ef. For Ef/T ^ 1, l|58|) 
can be simplified to 



E^f 
2T 



In L, 



(62) 



where we used (I61II is the last equality. Finally, the resistance of the array is given by i? = Rae^'-f with Uf = f{Ef, T). 
Expressing Ef as a, function of ly with the help of H62|l . we obtain the resistance of the short array in the form 



R{ly,T) 



Rq exp 



LT ^ \ agT 



1/2- 



agT 
21ni„ 



1/4 



agT 
21ni„ 



3/4 



(63) 



Expression H63I) determines the temperature dependence of the short array and is valid for the length of the array 
ly < e '^s^T ^ jn the opposite regime, ly > , the resistance is given by Eq. (|6()|l . 

E. Plasmon assisted resistance for short hopping time 

The integral over the time in H35|l should actually be cut off by the hopping time th. In the regime of very low 
temperatures and short hopping time over the break ^ tfi ^ the large argument expansion H49|) 

applies for m ^ while a short time expansion should be used in calculating Ij.{t). We obtain Ij.{t) = 1/2 + (-^^7^) 
from which it follows I±{uj) = ^S{uj). For Z{t), leaving only the factor with m = 1, we get 



Z{t) « 1 + 2k 



2LT 



COS —t + O 



(64) 



form which it follows 



Z{lu) w 27r(5(w) + 27rK 



2LT 



- / -KW \ J. / TTIW \ 



(65) 



Eq. H65|) shows that the bosonic spectrum remains discrete, corresponding to the spectrum of a single wire. Physically 
it means that the plasmons from the distant wires cannot reach the break during the hopping time, and the existence 
of plasmon bands plays no role for the electronic transport. In that regime, the influence of interactions reduces to the 
charging energy of a finite wire, which is the main assumption in the models of granular metals with high intergrain 
conduct anc o^^i^^ . Finally we note that the condition of short hopping time may contradict the assumption of an 
exponentially high resistance of a break. 



RESISTANCE AT HIGH TEMPERATURES 



At high temperatures, when T 3> the temperature broadening exceeds the interlevel separation in a single wire. 
To analyze the temperature dependence of resistance in that regime, it is convenient to rearrange the sum over m in 
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(|38|l with help of the Poisson summation formula 



CT— ±1 m— 1 



/ sinh[^ (m/3 + r)] 
I sinh (m/3 



da; In 



/ sinh (/3a; + r)] 



sinh 



2L 



{l3x + a)] 



- E 1- 

m— — oo 

2 / dx cos(27rfca;) In 




/ sinh + t)] 

l^sinh (/3a; + a)] 



(66) 



The last line of H66|l has been obtained using the Poisson summation formula. After rescaling the integration variable 
in the integral over x as y = ^^/3a;, we obtain 



S 



2L 

7r/3up J _ 



dy In 



/ sinh[y+^r] 
V sinh \y 



Integrating by parts in each term with a given fc, we finally obtain 



1+2Y: 



k=l 



(67) 



5 = 

fc=i 
1 
k 



cos 27rfc- 



/3 



cos 27rfc- 



coth 



27rfcL 



(68) 
(69) 



Let us rearrange the sum for S in form of a perturbative expansion around the result for a Luttinger liquid of infinite 
length. To this end we note, that if we replace coth = 1 for all k, then the sum over k reads 



k=l 



sin(7rT//3) 



sin(7ra//3) 



(70) 



which is the result for the bosonic correlator in the infinite length Luttinger liquid. Therefore, the expression for S 
can be written in the form 



S = log 



sin(7rT//3) 



sin(7ra//J) 



oc _^ 



k=l 



COS 27rfc 



cos 27rfc 



coth 



27rfcL 

/3up 



- 1 



(71) 



Here the first term is the infinite length result, and the sum over k denotes the corrections due to the finite length of 
the wire L. For large argument of coth, ^ ^^^^ ^ '3> 1, the /c-correction decays as j exp (^—^^^^ ■ Further calculation 

is performed for the zero order term in 1/L expansion. Evaluating the Fourier transform of Xi){t) 
performing the analytical continuation to real frequencies, we obtain 



-{KpS(T)) 



p and 



X,{c.) « -l^sin(W2)r(l - 



(72) 



The conductance of a single junction is then calculated by H36|l . Substituting w = —Ec in H72|l one can see that 
the linear current through the junction and hence its conductance does not fluctuate exponentially strongly with the 
charging energy Ec- Therefore, the resistance of the array is not determined by the resistance of a single break, it 
is rather given by the average over the resistances of all junctions. This result also suggests that at high enough 
temperatures the exponential suppression factor e'^"!'^ that is typical for the activated transport is compensated 
and even overridden by the exponential increase in the number of particle hole pairs that can provide the energy E^ 
to hop over the break. 

The fact that the leading contribution to the conductance at high temperatures is given by the result for infinitely 
long wires implies that the coherence of the single particle motion is broken already by a single hop between the two 
neighbor wires. The latter justifies the averaging over the all junctions in calculation of the resistance. In the result, 
the expression for the resistance can be written as a Drude formula a — e^NTfjm^ {N is the electron concentration 
and TO* is the effective mass) with the interaction and temperature dependent mean free time t/(T). The expression 
for the mean free time can be organized as an expansion in powers of the small parameter e ^ . 
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At a typical charging energy Ec ^ 2ttT, the derivative dX^/duj for u! = —Ec in (|36|l can be approximated as 



dXb 



_2a''sin(^) T{2~k)E, 



K-2 



(73) 



Thus, in that regime the leading term in the expression for the mean free time is temperature independent, Tf oc 
1/(E'^^'^). For comparison, phonon assisted transport in that temperature regime still has a thermally activated 
character with the preexponential factor T^^/^. 

At temperatures even larger than the typical charging energy, Ec < T, the expansion of 
parameter results in the following expression for the conductance of the single junction 



du 



in small 



r(i - k)v'{o)t''--^e, 



where V^'(O) is the derivative dV{x)/dx\x^Q of the hmction 

r(K/2 + a 



V{x) 



[7/>(k/2 + x) - V'(l - k/2 + x)] 



(74) 



(75) 



r(l-K/2 + a;) 

Subsequent averaging of the resistances of all junctions over the random charging energy E^ leads to the Drude 
expression for the resistance, where r/ exhibits the power-low temperature dependence typical for transport across a 
sliding Luttinger liquid, r/ oc T^~'^{E~^)^. Therefore, at high temperatures, the mean free time is determined by 
the Luttinger liquid interaction parameter n. The resistance for phonon assisted transport in that regime is given by 
a Drude formula with the logarithmic temperature dependence of the mean free time t/ oc InT. 



VI. PHONON ASSISTED TRANSPORT 



In this section we calculate the phonon assisted thermally activated transport in the system under consideration. 
For simplicity, we consider the case of noninteracting electrons, taking into account only the charging energy by the 
transfer of electrons between the wires. Then the single particle energy spectrum of each wire near the Fermi energy 
consists of equidistant energy levels with the interlevel distance ttvq/L, vq being the Fermi velocity. Furthermore, we 
assume the phonons to have a featureless density of states with the spectral density Z?ph = ojd / {2-kVs), where lou is 
the Debye frequency and Vg is the sound velocity. The phonon assisted transition rate between the wires and 1 that 
form the break can be calculated aa^ 

Toi = |ioi|'7o-ph ^ /(£„)[! -/(ez-ey)]ArB(i;e + e/-e„- el/), (76) 

where e„ = tivqu/L is the single particle energy with respect to the Fermi level, /(e) denotes the Fermi distribution 
and Nb{uj) denotes the Bose distribution, and 7e-ph is the electron-phonon coupling constant. The linear conductance 
of the break is given by 

CTbrcak = Tf^^Ol- (77) 



A. Low temperature regime 

Calculating the transition rate according to (|76|l we obtain 



r„i = |toi p7e-phe-''^= <! ^ + f^^^"' \, \ (78) 

1-e f^_^-p^\ 



Further calculation follows the procedure described in Section lllll and leads to the phonon assisted thermally activated 
resistance in the form 



-^phonon ~ a/ ^ e^^-^ . (79) 
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B. High temperature regime 

In the high temperature regime, T ^ T > Ec, we replace in l(7H)) /(cn) ~ 1 — /(e; — eV) « 1/2, and replace 
the summation over n, / by integration 

|ioip7c-ph f°° dx 
i 01 ~ - 



4 Jo e/^^^cg/'-L^^ - 1 

/-„, |2 /-.T^, 

ln|/3Sc|. (80) 



Attvq ' ' 47rt;o 



Substituting that result into H77|) . we obtain 

o-Qi w —\ii\PEc\. (81) 

As argued above, the resistance of the whole array is given by the average over the resistances of all junctions over 
the random charging energy Ec and the random hopping matrix elements tij , which results in 



-^phonon — 2 1 ^ l^ii l 1^ 



T 



(82) 



Therefore, at temperatures exceeding the typical addition energy, the resistance of phonon assisted transport has a 
logarithmic dependence of temperature. 



VII. CONCLUSION 



In this paper we calculated temperature dependence of resistance for plasmon assisted thermally activated transport 
in a particular model of a one-dimensional granular array. The considered model is relevant for experiments on 
transport across a highly anisotropic two-dimensional systems that can be represented as an array of parallel one- 
dimensional wiresiiiiSiiiii. 

We showed the existence of the regime, where the localization length of the plasmon excitations exceeds by far the 
single particle localization length. The latter allows plasmons to traverse the whole array and thus form an external 
bath of bosonic excitations necessary to suppress single particle Anderson localizationSi. Further peculiar property 
of the model lies in its closeness to the sliding Luttinger liquid model— as the description of interactions concerns. 
On one hand this property allows to use a powerfuU bosonization technique and to arrive at conclusions about the 
temperature dependence of the preexponential factor in the thermally activated resistance (|60|l . The calculation also 
gives an insight in the formation of multiplasmon complexes that create a bosonic bath with continuous density of 
states (|54f) . Moreover, the increase of the number of particle-hole excitation with temperature leads to the suppression 
of the thermally activated behavior at high temperatures, the resistance being given by the Drude formula with the 
mean free time that depends on temperature as a power law. The power low dependence of the Drude mean free 
time in the hight temperature regime is in general agreement with findings of the recent work*. At this point one 
might speculate about the analogy between that change of the character of transport and the recently reported finite 
temperature phase transition form the Anderson localized to a delocalized phase in interacting disordered system£:^^'«. 

On the other hand, our results show that the capacitive model used for the description of transport in granular 
arrays^ii^ii^ is inapplicable for that particular kind of array. This property again stems from the closeness of our model 
to the sliding Luttinger liquid one, where the relaxation time of charge density excitations is very long. The chosen 
specific model of a one dimensional array enabled us to show that thermally activated resistance has qualitatively 
different temperature dependence for plasmon assisted transport as compared to phonon assisted transport. 

Despite the specific quasi one-dimensional geometry of the grains in this model, the present investigation is believed 
to be of general importance for granular arrays with delocalized or weakly localized plasmons. 
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